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Abstract 

Frequencies in oscillating (3 Cephei stars are usually inferred by means of radial 
velocities measured from the Sim triplet AA 4552 - 4574 A. These lines, relatively 
insensitive to the variation of T e g through a pulsation cycle, show small equivalent 
width variations. 

In this study we aimed to verify if the behavior of radial velocities and equivalent 
widths measured from other ions are compatible with the one observed from Sim 
lines and than to verify the possible vertical stratification along the stellar atmo- 
sphere. For this reason we selected from our spectra a number of, unblended and 
well isolated, Cn, Nil and On lines besides the famous Sim triplet. All those lines 
cover the range in optical depth between -2.1 and -0.5. 

Unfortunately, we did not find any differences in the radial velocities behavior line- 
by-line and then we derived the frequency of the principal radial mode combining 
all the velocities derived from each spectral line separately. The inferred frequency 
was fi = 5.249677 ±0.000007 c d _1 . 

Another important task we would like to accomplish with this paper is to make 
available to the community our large sample of spectroscopic data, that is 932 
velocities and equivalent widths measured from our sample of C, N, O and Si lines. 
All the spectra were acquired at the 1-meter class telescope of the stellar station of 
the INAF - Osservatorio Astrofisico di Catania, in the period starting from July, 
the 27 th 2005 to November, the 1 st 2006. 
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1 Introduction 



The star (3 Cephei (HD 205021) is the prototype of a class of hot pulsating 
variables. For complet e reviews on (3 Cephei stars, we refer the readers to 
Aerts fc De Cat (2003), wh o desc ribe in detail the line-profile variations and 



to ISterken fc Jerzykiewica (119921 ). who give an overview of the photometric 
behavior of all the j3 Cep stars known up to that time. 

From the spectroscopic point of view, the data used to recover frequencies 
of oscillations are commonly the radial velocities. Checking through the most 
recent literature published on this topic, it is easy to realize that the common 
guideline is to derive the velocities from the lines of the Sim triplet at 4552, 
4567 and 4574 A. There are three common reasons for this choice: they are 
strong lines not affected very much by blending and almost insensitive to the 
temperature variations (IDe Ridder et all |2002| ). 



Up to the study by lAerts et al.l (11994 ). (3 Cephei has been believed to os- 
cillate with only a single radial mode. In their analysis these authors re- 
vealed the multi-periodicity of this star detecting three frequencies in to- 
t al: fi =5.24 9 7104 c d" 1 , f 2 = 5.385 c d" 1 and f 3 = 4.920 c d" 1 . Later on, 



Telting et al.l (119971 ) by a CLEANing analysis of the same data added two 
more frequencies: / 4 = 5.083 c d _1 and f 5 = 5.417 c d -1 . 

Further, the star j3 Cep (V = 3.2) is actually a complicated multiple system. In 
fact the star, beside to have a visual companion (V = 7.9) at a distance of 13.4" , 
is also a membe r of a spectroscopic system whose second star was discovered 
by iGezari et al.l (Il972l ). using speckle interferometry, at a distance of ^0.25" 
(V = 6.6). The parameters of t he close binary orbit have been determined later 
by iPigulsky fc Boratynl (119921 ) from speckle interferometry and the variations 
in the pulsation period, due to the so- called light time effect. Recent speckle 
measurements by lHartkopf et al.l (120011 ) place the position of the companion to 
about 0.1" from the primary. More observation are still necessary to improve 
the orbital solution. 



Another question, people are still debating, is the nature of the random vari- 
able H Q emission which was the cause of the Be status assigned to (3 Cep. 
The problem was that typically Be stars are rapid ro t ators while (3 Cep has 
an equ atorial velocity of 26 km s _1 ( iMorel et al.l . 120061 ) . lHadrava fc Harmanec 
( 119961 ) . attempted to solve this problem argued that the observed emission is 
due to the secondary component of the spectroscopic binary which is proba- 
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bly a Be star. They, for the first time, se parated the emission and absorption 



component of the H a profile. Recently, ISchnerr et al.l (120061 ) confirmed this 
result using data obtained at Nordic Optical Telescope. Unfortunately, they 
obtained only one spectrum and, then, no conclusion about the variability of 
the emission could be drawn. 



Starting from the pioneering work by lOsakil (Il97ll ). a lot of studies have been 



undertaken with the aims to infer the oscillation modes from the analysis of 
line profile variation (LPV). It is generally accepted that the weakness of these 
studies is the lack in the knowledge of how some important thermodynamic 
quantities change during the pulsation. For instance, during each pulsation 
cycle effective temperature, surface gravity and microturbulent velocity are 
expected to change with time and than the knowledge of their effect on the 
line profiles could avoid errors in the modes determination. 

It is a matter of fact that, because of its low rotational velocity, the spectrum 
of (3 Cephei is full of lines generated by various chemical elements at different 
depths in the photosphere. The analysis of the radial velocity inferred from 
a sample of lines opportunely selected could be useful to study the velocity 
stratification along the photosphere. A similar approach has been tempted by 



Baldry et al.l (119981 ) in a spectroscopy study of the roAp star a Cir. In this 
study latter authors selected a sample of spectral region containing various 
metallic lines and they found a dependence of the amplitudes with the depth 
along the atmosphere. 

With the aim to infer if a similar stratification is present also in (3 Cep, we 
selected from our spectra as many as possible lines not affected by blends, 
sufficiently isolated to allow an easy measurements and well distributed along 
the atmosphere. A complete list of selected lines is reported in Tab. [2J In this 
paper we present a large time-resolved dataset of line-by-line radial velocities 
collected from July 2005 to November 2006, for a total of 932 spectra and, 
moreover, for the first time we present also their equivalent widths. 

Velocities and equivalent widths can be retrieved from the CDS via anonymous 
ftp. 



2 Observations and data reduction 



All the data analyzed in this study have been obtained at the 91 cm telescope 
of the INAF - Osservatorio Astrofisico di Catania. The telescope is fiber linked 
to a REOSC echelle spectrograph, which allows to obtain R = 20 000 spectra in 
the range 4300-6800 A. The resolving power has been checked using emission 
lines of the Th-Ar calibration lamp. Spectra were recorded on a thinned, back- 
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Table 1 

Journal of observations. JD are in the form JD - 2450000, n is the number of spectra 
collected in each night. 



JD 




JD 




3579 


56 


3933 


49 


3580 


44 


4013 


35 


3581 


46 


4016 


30 


3635 


65 


4017 


26 


3672 


64 


4018 


41 


3673 


95 


4019 


43 


3674 


82 


4021 


19 


3871 


58 


4022 


18 


3872 


61 


4023 


5 


3873 


64 


4041 


31 



illuminated (SITE) CCD with 1024 x 1024 pixels of 24 fim size, typical readout 
noise of 6.5 e~ and gain of 2.5 ph/ADU. 

Our observations (932 spectra in total) were spread over a baseline greater than 
one year, precisely from July, the 27 th 2005 ( JD = 2453579) to November, the 
1 st 2006 (JD = 2454041), for a total of 462.14 days. The exposure time was 
set to 5 minutes, with the exception of a couple of nights in October 2006 
during which we had to increase it to 10 minutes because of the bad weather 
condition. This leads to a temporal resolution of ~ 2% and ~ 4%, respectively, 
of the pulsational period, small enough in both cases to avoid phase smearing 
effects. The signal-to-noise ratio of our spectra resulted always above 150. 

The stellar spectra, calibrated in wavelength and with the continuum nor- 
malized to a unity level, were obtained using standard data reduction proce- 
dures for spectroscopic observations within the NOAO/IRAF package, that 
is: bias frame subtraction, trimming, scattered light correction, flat-fielding, 
fitting traces and orders extraction and, finally, wavelength calibration. IRAF 
package rvcorrect has been used to include the velocity correction due to the 
Earth's motions, all the spectra were then reduced into the heliocentric rest 
of frame. 



Systematic errors on radial velocities have been estimated observing stars 
with constant and well known r adial velocity taken from the list of standard 
stars published bv lUdrv et all fjl999h : HD 12929, HD 20902, HD 186791 and 
HD 206778. Whenever it was not possible to obtain a spectrum of a standard 
star, we evaluated the systematic shift from the wavelength positions of the 
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Fig. 1. Phase diagrams of the radial velocity variations derived by carbon, nitrogen, 
oxygen and silicon lines (left side) and velocities computed from average over all 
the lines selected (right side). Parameters of the over-imposed sinusoidal fits are 
reported in Tab. [2j Error bars extend for ± 1 a. 



strong telluric lines between 6275 A and 6320 A. 



3 Radial velocities and equivalent widths measurements 



As we stated before, radial velocities in this class of pulsating stars are usually 
measured from the Sim triplet. Here we have considered the heliocentric V ra( j 
measured from lines of other ions with the final aim to verify if their behavior is 
different line-by-line. Thus, we selected all the strong, unblended lines present 
in the spectral range covered by our spectra. These lines are reported in the 
first column of Tab. [2] and are, precisely, 2 lines belonging to Cll, 5 to Nil, 2 
to Oil and 3 lines to Sim. 

For each of these lines we adopted as radial velocity the value obtained from 
the barycentric wavelength which has been converted in radial velocity via the 
classical Doppler shift formula. The integrations evaluated to compute Ab ar 
and the equivalent width (EW) of each line extend over a carefully chosen, 
constant region in wavelength and they have been evaluated after a local re- 
normalization of the continuum close to the edges of the line. This procedure 
let us to minimize the systematic effects due to the arbitrary choice of the 
lines boundaries or the position of the continuum. Errors on V ra d and EW 
have been computed from a formal application of the error propagation rules 
to the used equations, assuming that the main source of error is the noise in 
the observed spectrum. 



5 



Table 2 

Parameters and relative errors of the sinusoidal fits presented in Fig. Q] and in 
Fig. [3j For each line we report respectively: wavelenght, log To, derived frequecy of 
the principal radial mode (in parentheses the estimated error), 70, semi-amplitude 
and phase of the variation. Columns from seventh to tenth reported the parameters 
derived for the fits of EWs, that is: relative amplitude, average, semi-amplitude and 
phase of the variations. 







HRV fit 


parameters 




EW fit parameters 


A 


log r 


frcq 


70 


K 


<t> 


rcl. 


70 


K 


<t> 


(A) 




(c/d) 


(km s _1 ) 


(km s _1 ) 




(%) 


(mA) 


(mA) 




Cn 6578.052 


-2.11 


5.24967(1) 


-3.2 ±0.6 


15.5±0.1 


0.526 ±0.001 


9.5 


88.2 ±0.2 


8.4±0.2 


0.28 ±0.01 


Cn 6582.882 


-1.84 


5.24967(1) 


-3.3±0.5 


13.5±0.1 


0.544 ±0.002 


9.9 


82.4 ± 0.2 


8.1 ±0.2 


0.19 ±0.01 


Nil 4607. 153 


-0.53 


5.24969(2) 


-1.9±0.4 


15.3±0.2 


0.456 ±0.002 


9.1 


56.0 ±0.2 


5.2±0.3 


0.19 ±0.01 


Nil 5010.621 


-0.53 


5.24967(2) 


-1.8±0.5 


13.1 ±0.2 


0.530 ±0.002 


7.4 


51.0 ±0.3 


3.8±0.4 


0.24 ± 0.02 


Nil 5666.629 


-0.97 


5.24967(2) 


-3.1±0.5 


14.1 ±0.2 


0.528 ±0.002 


8.4 


83.3 ±0.2 


7.0±0.2 


0.18 ± 0.01 


Nil 5676.017 


-0.82 


5.24967(2) 


-4.3±0.4 


13.3±0.2 


0.536 ±0.003 


7.9 


69.4 ±0.2 


5.5±0.2 


0.19±0.01 


Nil 5679.558 


-1.27 


5.24968(1) 


-4.5 ±0.6 


14.3 ±0.2 


0.500 ±0.002 


7.1 


120.2 ± 0.2 


8.5±0.2 


0.21 ±0.01 


On 4591.010 


-0.82 


5.24967(1) 


-5.3±0.7 


15.1 ±0.1 


0.544 ±0.001 


1.2 


120.5 ±0.3 


1.5±0.4 


0.05 ±0.04 


On 4661.643 


-0.97 


5.24966(1) 


-4.7±0.7 


15.9±0.1 


0.597 ±0.001 


1.9 


121.5 ±0.3 


2.4±0.4 


0.39 ±0.03 


Sim 4552.622 


-1.70 


5.24967(1) 


-1.9±0.8 


14.4 ±0.1 


0.535 ±0.002 


2.5 


209.0 ±0.2 


5.2±0.3 


0.33 ±0.01 


Sim 4567. 840 


-1.41 


5.24967(1) 


-2.6±0.7 


15.0±0.1 


0.556 ±0.001 


2.7 


174.3 ±0.2 


4.8±0.3 


0.36 ±0.01 


Sim 4574.757 


-0.97 


5.24968(1) 


-3.2±0.6 


15.4±0.1 


0.505 ±0.001 


2.9 


117.9 ± 0.2 


3.4±0.3 


0.27 ± 0.02 


All 


5.249677(7) 


-3.5 ±0.4 


15.5±0.1 


0.516 ±0.001 





Single values of radial velocities and equivalent widths are presented only in 
the online tables^, available at the CDS, which are organized as follow: the 
first column reports the heliocentric julian date, the others reported line-by- 
line the measured radial velocities and equivalent widths together with their 
errors. Blank space means a discarded line profile contaminated by cosmic 
rays. 



4 Data Analisys 



One of the most used code to analyse a large set of periodic data in order to 
obtain frequency, velocity of center of mass (70), amplitude of the variation 
(K) and phase ((f)), is Period04 developed by lLenz fc Breguerl (120051 ). We used 
this code to perform a line-by-line sinusoidal fit of our velocities. For each 
selected line, we reported in Tab. [2] frequency, 70, amplitude K and phase 
cf) and relative errors as determinated by using that code. Each box on the 
left side of Fig. [1] shows the velocities derived line by line folded with the 
corresponding frequency and with the sinusoidal fit over-imposed. 



With the aim to search for vertical stratification, we plotted all the quantities 
reported in Tab. [2] regarding the velocity vers us optical depths. To com pute the 



optical depths, we used the code XLINOP (IKurucz fc Avrettl . Il98ll ) applied 



1 Table 1 reports lines from Cn 6578 to Nil 5676, Table 2 reports lines from Nil 
5679 to Sim 4574 
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5.24970 
5.24968 
5.24966 
5.24964 



Fig. 2. In this figure we report as a function of optical depth the computed frequen- 
cies, center of mass (70), amplitudes and phases derived from the fit of velocities for 
each of the selected lines. Meaning of the symbols is: circles (red) carbon lines, stars 
(magenta) silicon lines, triangles (blues) oxygen lines and boxes (green) nitrogen 
lines. 

to LTE atmospheric model with T e ff = 24000 K, log g = 3,9, solar metalicity 
and £ = 8 km . These atmospheric values have been derived by us in the 
framework of a detailed abundances revision that it is currently under study. 
As we shown in Fig. [2j there is no evidence of some kind of stratification 
through the stellar atmosphere, at least in the range of optical depth probed 
by our lines. 

Since the results obtained line- by-line are consistent, with the aim to improve 
the quality of our measurements, we computed the radial velocity as an average 
over the values obtained from all lines, without paying any attention to the dif- 
ferent chemical species. In this case, the fitting procedure described before gave 
us the results presented in the last row of Tab. [2] and then we adopted it as fre- 
quencies of the radial mode of (3 Cep: f\ = 5.249677 ± 0.000007 c d _1 . The cor- 
responding phase diagram is shown in the right side of Fig. [U The epoch used 
for the plot is that of the minimum radius, that is HJD( Rt min ) = 24 53579.2649. 



Since the study by lPigulsky &: Boratynl (119921 ) . it is well known that (3 Cephei 



shows apparent changes of its pulsation period because the light -time effect 
induced by the orbital motion of the binary system. In our conclusion we have 
neglected this effect since our time baseline covers scarcely 1.4 % of the orbital 
period, estimated by those authors in 91.6 yrs. 
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Nil A5679.558 A 




Fig. 3. (Left side) Variability of the measured equivalent widths for all the selected 
lines. Parameters of the sinusoidal fits over-imposed are those reported in Tab. [2J 
(Right side) Equivalent width variations for the Nil A5679 A observed in ten different 
nights. In each box we indicated the julian day - 24 50000. 



As a by-product, we measured also the equivalent widths of each selected 
line. In Fig. [3] (left side), we show the phase diagrams obtained for all the EW 
measured for the selected lines phased with the ephemeris derived from the ve- 
locities. As we stated in the previous section, all these lines are unblended and 
sufficiently isolated to avoid contamination from nearby lines. Sinusoidal fits 
over-imposed on the experimental data have been computed neglecting those 
points for which their relative errors exceed 20 % and keeping all the points 
within ± 3 a from the mean value, where a is the standard deviation. The pa- 
rameters of the sinusoidal fits with their errors are presented in Tab. [2j We did 
find variations with small amplitude for some lines as for example Oil 4661 A 
and Sim lines (relative amplitudes <3%). A more marked variability, up to 
~ 10 %, has instead been inferred for Nil and Cll lines. 



Despite difficulties in measuring the EW of a spectral line, for instance the 
choice of line limits and continuum, the phase diagrams showed a scatter not 
justified if related to errors. For the sake of clarity, we showed in Fig. [3] (right 
side) the EW measured in 10 different nights for the Nil A5679 A line. Similar 
behavior have been observed also in the variations of the others spectral lines 
here considered. A so complex changing, obs erved night by nigh t , could be as- 



cribed to the multi-periodic nature of (5 Cep (lAerts et ajj (119941 ). iTelting et al. 
fll997h l 
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5 Conclusions and discussion 



Spectroscopic observations of (3 Cephei have been carried out at the stellar 
station of the INAF- Osservatorio Astrofisico di Catania, in total we collected 
932 spectra with a time baseline covering more than one year. 

As for the radial velocities we checked the behavior of velocities coming from 
different lines and for each group of velocities we computed the pulsation 
frequency of the radial mode, the 70, amplitude and phase of the relative 
sinusoidal fit. Plotting these values as a function of the optical depth of each 
line we do not see any evidence of stratification through the stellar atmosphere. 
For what that concerns equivalent widths, we showed their behavior in time 
and, in the case of the Nil A 5679.558 A, the possible influence of the non- 
radial pulsation on the amplitude and shape observed night by night. Further, 
we did find clear variability for Cn and Nil lines, while Sim and On A 4661 A 
show very small amplitude. For On A 4591 A we concluded that its EW is 
constant with time, at least at the resolving power of our spectra. 

At the end, we think that such large data set of data, radial velocities and 
equivalent widths, could be useful for a number of studies, like for instance 
the study of the orbit utilizing the light-time effect or frequencies analysis, so 
we would like put them at astronomical communinity's disposal. 
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